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ABSTRACT 

Ha equivalent widths and UBV magnitudes are calculated for Be star disk models that grow in size 
and/or density with time. We show that these simple models are consistent with the known Be star 
classes of positive and inverse correlations between long-term variations in Ha and V magnitude as 
identified by Harmanec. We support the conclusion of Harmanec that the distinction is controlled by 
the inclination of the disk to the line of sight. We demonstrate that the strength of these correlations, 
particularly those of an inverse correlation where the system becomes fainter as the Ha emission 
strength grows, is strongly influenced by the scale height of the inner Be star disk and the extent of 
the gravitational darkening of the central B star due to rapid rotation. This dependence may allow 
coordinated spectroscopic and photometric observations to better constrain these poorly known Be 
star parameters. 

Subject headings: stars: circumstellar matter - stars: emission line, Be 



1. INTRODUCTION 

Be stars are defined observationally as non-supergiant 
B stars that have shown emission at le a st on ce in the 
Balmer seri es of hydrogen Jj aschck ct at (198l]h s ee als o 
reviews bv iKogure fc Hiratal 1(19821) : iSlettebaH ((1981 ): 
iPorter fc Riviniusl (J2003[ r)~ Be stars also exhibit an 
IR e xcess relative to that expected from a normal B 
star (|Gehrz et al.lll974t iCote fc Waters! 1 1 98 7D and linear 
continuum polarization of up to 2% ( McLean fc Browril 
119781: iPoeckert et al.l 119791 ) . All of these observations 
are consistent with the presence of an equatorial disk 
of gas surrounding the central B st ar, an idea ("a nebu- 
lous ring of gas") first suggested bv iStruve] ((1931). This 
basic picture now has ample direct interferometric evi- 
dence, and the physical model of a flattened, equatorial 
disk surrounding the centra l B star is well esta b lished 
([Dougherty fc Tavlorl 119921: lOuirrenback et all 119931: 



Stee et all 119951: lOuirrenback et al.l 119971: iTvcner et al l 
20051: iSteel 120111) . About 17% of all non-supergiant B 



stars in the Milky Way are Be stars, although the frac- 
tion varies widely with spectral type, reach ing a maxi- 
mum of ~ 34% at Bl dZorec fc Briotl 119971 ). Recently, 
iTarasov fc Malchenkol (|2012l ) have found that the frac- 
tion of early-type Be stars (B0-B3) peaks in Galaxy clus- 
ters of ages between 12-20 Myr, suggesting that evolu- 
tionary status plays an important role in the Be phe- 
nomena. It also seems well established that Be stars are 
more frequent in lower metallicity environments, such 
as the Magellanic Clouds (IW isnicwski & Biorkmarj |2006l : 
IMartavan etaIll20M l2007bL 120101) . 

Be stars exhibit variability on a wid e variety of time 
scales, from minutes to years /decades (jHarmaned 1 1 9831 : 
IPorter fc Riyirmisl [200l . The origin of short term vari- 
ability is usually attributed to pulsations within the B 
star photosphere (|Baadd 12055 iHuat et all [20091. and 
there is speculation tha t pulsations may play a role in the 
disk formation (jOsakil Il999t !Cranmeill2009l ). However, 
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corotating structures, such as photospheric "clouds," 
have also been p roposed to e xplain these variations 
(Harmanec] [l999t iBalonal [25551) . Lon ger term varia- 
tions seem connected with the circumstellar disk, either 
the formation and dissipation of the disk itself, or hy- 
drodynamical processes occurring within the disk, such 
as spiral d ensity waves which lead to the cyclic V/R 
variations (lOkazakil Il997t iHummel fc Hanuschikl 119971: 
iCarciofi et al.ll2009D . 

One of the striking aspects of Be stars is episodic cre- 
ation and destruction of the disk for some Be stars. This 
is exhibited as a transition in the observed spectrum from 
a Be star to a normal B star and/or to a Be shell star. 
Such events present a unique opportunity to study both 
the structure of the disk and the physical processes re- 
sponsible for its creation. 

Even though the process of disk formation is not well 
understood, there are observations t hat seem to show 
disk formation and dissipation (eg. see iWisniewski et al.l 
120101 ). Many Be stars show outbursts, characterized by 
rapid variability in the emission line profiles. Such vari- 
ations are normally i nterpreted as events of rapid, dis- 
crete mass ejection (jRivinius et al.l fl998al Tbl). although 
a definitive observation of material escaping the stellar 
surface seems lacking (Harmanec, 2012, private commu- 
nication). These outbursts are thought to be the re- 
sult of material that is ejected from the star and are 
mostly seen in early-type Be stars. Stars such as uj CMa 
(B2IVe), FV CMa (B2IV-Vne), FW CMa (B2Vne), 
H Cen (B2Vnpe), ry Cen (B1.5Vne) and 28 Cyg (B2.5Ve) 
are examples of early t ype Be stars with outburst activity 
(|Rivinius et al.ll2001al) . 

Some Be stars are seen to repeat outbursts regularly. 
In these cases, the disk is replenished with material ev- 
ery few months to few years, and the outbursts fill up 
the inner region of the disk with the material. If an out- 
burst is not followed by another one soon, the material 
falling back on the star, or propagating outwards in the 
disk, will create a low density region near the star. The 
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longer the time before the next outburst, the emptier 
the region becomes. This empty region between the star 
and the inner edge o f the disk makes the disk structure 
look more like a ring ((Mcilland ct al. 2006). This "inside- 
out" mass loss can be seen in many s tars, for example 
7] Cen (jRivinius et al.l[200Tal ). 60 Cyg (jWisniewski et al.1 
I2010D . and 28 CMa (jStefl et al.lll995l ). where observers 
measure the variability of the width at the base of the 
e mission lines a nd / or spectropolarimetric data. 

IClark et al.1 (|200l give observations for o And 
(B6IIIpe) which show the formation and destruction of 
the disk in an inside-out manner. In their 15 year data 
set, the star goes from a normal B star to a Be shell 
star and back to normal B star, with the cycle repeating. 
This transition takes around 700 days. Short term out- 
bursts have be en noticed which r eplenish the disk with 
more material. IClark et al.l ([2003D model one of the shell 
episodes, and it is concluded that there is a decrease in 
density in the innermost regions of the disk. This de- 
crease is shown to slowly propagate outward, indicating 
an inside-out mass loss in the disk. 

Why some B stars, and not others, become Be stars 
is currently unclear. Rapid r otation of the central B 
star seems to play a key role (porter fc Riviniusl [2003), 
but details are still lacking as a definitive determination 
of the actual rotation rates of the central stars is com- 
plic ated by the potential effects of gravitational darken- 
ing (Townsend et al.ll2004l: iCranmen 120051 : iFremat et al.l 
2005). The evolutionary status of the central B star 
may also play an importan t role (jEkstrom et al.l [20081 : 
ITarasov fc Malchenkol [20121 ) . 

Also unclear is the exact mechanism(s) that creates Be 
star disks. Keplerian rotation, now established for Be 
star disks ([Hummel fc Vranckenl 120001 : lOudmaiier et all 
2008), sugg ests viscous decretion as the correct physi- 
cal process (|Lee et al.lll991t IPorterl lT999[). but how ma- 
terial is feed into the inner edge of the disk is un- 
known (jCranmerl 120091 : IKrticka et ail 1201 ID. It impos- 
sible pulsation plays the critical role (jCranmerl [2009). 
and there is ta ntalizing observationa l evidence support - 
ing this view (jRivinius et all 119981: iHuat et al.l I2009D ■ 
The rol e of binarity in the Be phenomena is also 
uncle ar ([Porter fc Riviniusl 120031 : ITarasov fe~M alchcnko 
12012ft . although the majority of Be stars do not seem 
to be the result of binary evolution ([Baadd 119921 : 
iVan Bever fc Vanbeverenlll997D . 

2. OBSERVED Ha AND UBV CORRELATIONS 

IHarmanecl (|1983l) describes two classes of long-term 
variation for Be stars that can be established from con- 
temporaneous observations of photometry (such as visual 
magnitude and UBV colours) and a measure of the Ha 
emission strength (such as its equivalent width or peak 
continuum contrast): positive correlations and negative 
or "inverse" correlations. The distinction between posi- 
tive and inverse corr elations is thought to be geometric 
(see IHarmanecl Il983[ ) . as will be described shortly. 

A positive correlation is seen for most Be stars for 
which the relevant observations are available, and it is 
characterized by a decrease in the (visual) magnitude of 
the system with increasing strength of the B aimer Ha 
emission (such as an increase in its equivalent width in 
emission). In a (U — B) vs. (B — V) colour-colour dia- 
gram, Be stars with a positive correlation change their 



lumin osity class but not their spectral type (Harmancc 
1983). Increasing Ha signifies an increasing disk, and it 
therefore seems natural to expect an increase in the over- 
all brightness of the system, star-plus-disk, as the disk is 
built. 

Many stars are kno wn to exhibit th i s positive correla - 
tion: k Dra (B6 IIIpe:lSaad et alj|2004t Uuza et aT] |1994D, 
60 Cyg (BlVe: IKoubskv et al.l 120001) OT Gem (B2Ve ; 
IBozic et al.lll999D , V442 And (B2IVe: IBozic et al.l[200l , 
EW Lac (B3IVe: iFloquet et al.l I2Q00D tt Aqr (BlVe; 
Nordh fc Qlofssonlll977l). QR Vul (B3Ve: I Pavlovski et al.1 
1983h. a Cen (B2 Vnpe; IPachs et allll992h. 6 Per (B2Ve: 
Bozicet al.l [l995h . -y Cas (BOIVe: iDoazon et al.l Jl9"83|h 



28 Tau (B8IVev: Iffirata et aLlfgOOOfc iTanaka et alj|2007li . 
and uj CMa (B2IVe: IHarmanedll998h . ~ 

On the othe r hand , an inverse correlation, as described 
by [Harmancc (1983(), is characterized by an an increase 
in the magnitude of the system with an increase in the 
strength of the Balmer emission. In the (U — B) vs. {B — 
V) colour-colour diagram, Be stars that show an inverse 
correlation move along the main sequen ce, changing thei r 
spectral type, but not luminosity class ([Harmanedl 1 983h . 

Fewer stars are known to exhibit th is inverse correla- 
tion: 4 Her (B9e: IKoubskv et al.l ll997h. 88 H er (B7pshe; 
IDoazan et all [1982J), V129 4 Aql (BOVe; |H ^rn_e£aL| 
EMI, and T] Cen (B1.5IVne: IStefl et al.l[l99^h . 

An inverse correlation is thought to occur when a Be 
star disk is viewed more ed ge-on (i.e. at higher inclination 
angle i) (jHarmaned 1 1 983f ) . The forming disk can then 
act to reduce the brightness of the system by blocking 
the light from the stellar disk, while the small projected 
area of the disk on the sky keeps the disk emission to 
a minimum. If the critical inclination angle required to 
observe an inverse correlation is large enough, inverse 
correlations will be statistically less likely to be found 
than positive correlations. 

Hence the basic picture is that for Be stars observed at 
i < icrit, disk formation should result in a positive corre- 
lation, while Be stars observed at i > z cr i t should exhibit 
an inverse correlation. Of course, i cr a may depend on 
spectral type and/or the particulars of how the disk is 
being built. However, i cr n must depend to some extent 
on the thickness of the disk. Hence relative numbers 
of positive and inverse correlations may be able to con- 
strain, in a general sense, the thickness of Be star disks. 
In fac t this argument was used originally by Ha rmanecl 
(1983) who noted that the occurrence of positive and 
inverse correlations in the long-term variability was, at 
the time, among the strongest arguments for flattened 
Be star disks. This conclusion is strengthened by con- 
sistency of the correlation classification: if the record of 
observations for a given Be star is long enough, the same 
type of correlatio n (positive or inverse) will be repeated 
(jHarm ancc 200(3) ■ As the probability of observing a Be 
star with inclination i is just sinidi, the probability of 
observing an inclination i > z cr i t is cos z cr it . Using the 
number of cited positive and inverse correlations above, 
4 out of 16 Be star systems show an inverse correlation 
suggesting i cr it ~ 75°. However, this is only suggestive 
because the stars used to estimate the fraction were not 
selected from a homogeneous sample but simply culled 
from the literature. 

To date there have been no systematic calculations 
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demonstrating these correlations, either positive or neg- 
ative, as a function of viewing angle i. Of particular 
interest is the case of inverse correlation. At high in- 
clination angles (i near 90°), the disk is viewed edge-on 
and the blockage of starlight is caused by the thickness 
of the disk. While almost all models of Be star disks as- 
sume very thin disks with scale heights set by hydrostatic 
equilibrium, there are few direct constraints of this as- 
sumption, and several line s of evidence suggest that the 
disks m ay be much thicker (|Arias etai]2 006: Z orec et all 
2007b). In addition, the geometry of the inverse correla- 
tion, namely blockage of the equatorial regions of the 
star, suggests that any gravitational darkening of the 
stellar surface may also play a critical role. Hence the 
interplay of disk scale height and gravitational darkening 
may allow new constraints to be placed on these param- 
eters by the detailed study of inverse correlations. It is 
the purpose of this paper to provide a preliminary in- 
vestigation of these correlations using radiative transfer 
models for Be star disks. 

3. CALCULATIONS 

The thermal structure of the model Be star disks used 
in this work was co mputed with the bedisk code of 
ISigut fc Jones! (|2007D . This code enforces radiative equi- 
librium in a photoionized, circumstellar disk of a pre- 
scribed density structure by including heating and cool- 
ing processes for the nine abundant elements (H, He, 
CNO, Mg, Si, Ca, & Fe), each over several ionization 
stages. Details of the atomic models and atomic data, 
as well as an overvie w of the bedisk code, are given in 
ISigut fc Jonesl IpOOl . The main energy input into the 
disk is assumed to be the photoiqnizing radiation from 
the central star. Th e older lKuruczl (fl99 3) LTE photoion- 
izing fluxes used by ISigut fc Jonesl (120071) were replaced 
with t he newer, non-LTE calculations of lLanz fe Hubenvl 
(2007). To compute the Ha equivalent widths and UBV 
magni tudes of the models, the BERAY code of iSigutl 
was used. In all cases, this modelling approach 
considers only an isolated Be star and neglects any po- 
tential effects of a binary companion which are outside 
the scope of the current work. 

The density structure of the disk as a function of time 
was assumed to be of the form 

p{R,Z,t)= Po {t){^y e-^)\ (1) 

where R and Z are the cylindrical co-ordinates for the 
axisymmetric disk, and i?* is the radius of the central 
B star. The quantities n and p (t) are adjustable pa- 
rameters that fix the density structure of the disk. The 
disk was assumed to extend from the stellar photosphere 
at R = i?„ at the inner edge to an outer disk radius of 
R = Rd(t). This simple density model has been very 
successf ul in interpreting a wide range of Be star obser- 
vations (jGies et al.1120071 : iTvcner et al.H2008t Uones et al.1 
2008). The disk is assumed to be in Keplerian rotation 
and hence rotationally supported in the R direction. 

The vertical (or Z) dependence of Eq. ([1]) contains the 
scale height function H which is defined as 

( 2R*kT HE \ 1/2 
\ GM» fimnj ' 



Here M* is the mass of the central B star, p is the 
mean-molecular weight of the disk gas, and The is an as- 
sumed isothermal disk temperature. The Z-dependence 
of Eq. (JT]) and Eq. ([2]) follow from the assumption that 
the disk is in vertical hydrostatic equilibrium set by 
the Z component of the stellar gravitational acceleration 
and an assumed isothermal disk temperature, The- In 
Eq. ([2]), the mean molecular weight, /i, was taken to be 
0.68, appropriate for a fully-ionized hydrogen gas with 
a 10% mixture of neutral helium. The isothermal disk 
temperature, The (which is taken to be a constant for 
the entire disk) is further discussed below. 
Eq. can be instructively written as 



where c s is the local sound speed in the disk (set by tem- 
perature The) and Vk is the Keplerian orbital velocity 
at R. As the radiative equilibrium temperatures in the 
disk give sound speeds on the order of ten kms -1 while 
the orbital speed is many hundreds of kms -1 , the hydro- 
static model predicts disks that are very thin, often only 
few percent of the stellar radius near the star. 

Note that The is only used in Eq. © to fix the ver- 
tical density structure of the disk; the actual tempera- 
tures in the disk, T(R, Z), are found by enforcing radia- 
tive equilibrium. ISigut et all (|2009D consider consistent 
models in which the vertical structure of the disk is found 
by integrating the equation of hydrostatic equilibrium in 
a manner consistent with the radiative equilibrium disk 
temperatures: this treatmen t eliminates the ne ed for the 
parameter The- In general. fSigut et ail {2009) find that 
vertical disk scale heights are generally overestimated by 
Eq. @, using the typical value of The ~ 0.6T e ff, in the 
inner regions of the disk due to a cool zone that tends 
to form in the equatorial plane of the disk near the star 
for disks of sufficient density (i.e. po). Nevertheless, as 
will be shown in a later section, disks using Eq. ([2]) with 
The = 0.6 T c ff are already sufficiently thin to have sig- 
nificant consequences for the strength of the predicted 
Ha-V magnitude correlations Q 

The thermal disk models are thus described by the 
spectral type of the central B star (which is assumed to 
fix the stellar mass, radius and effective temperature) 
and the parameters in Eq. ([TJ that fix the density of the 
disk: p , the base density (in gcm~ 3 ), n, the power-law 
index, and i?d, the outer disk radius. The isothermal disk 
temperature used in Eq. ([2]), The, is also a parameter 
of the models, and variations of this parameter are the 
subject of Section l4~2l 

In the present work, two spectral types were consid- 
ered for the central B stars: B1V and B5V, representing 
an early and later-type Be star. The assumed funda- 
mental parameters are given in Table [TJ Be stars are, 
of course, known as rapid rotators and hence the role of 
gravitational darkening of the stellar photosphere must 
be considered. This effect will be explicitly considered in 
Section [O] 

1 IMc Gill ct al. (2012) demonstrate that the effect of gravita- 
tional darkening on the disk tempera tures can result in scale 
heights even smaller than those found by Sigut ct al. (2009). Hence 
the use of Eq. ([2} in the current work should be considered an upper 
limit to the thicknesses of disks in hydrostatic equilibrium. 
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Table 1 

Parameters adopted for the central B 
stars. 



Parameter B1V B5V 

Mass (M ) 13.2 5.9 

Radius (Rq) 6.4 3.9 

Luminosity (L ) 1.4 ■ 10 4 6.9 • 10 2 

T off (K) 25 000 15 000 

log(g) 4.0 4.0 

The 15 000 9 000 



Note. — The mass and radius cali- 
brations are taken from Cox (2000). 

As this paper is concerned with the optical signature 
of disk formation and dissipation, some prescription for 
these processes is required. One-dimensional, hydrody- 
namic simulations that evolve the disk surface density 
as a function of time within the framework of viscous 
decretion have reached a high deg ree of sophistication 
(Carciofi 20101 : lHaubois et al.l I2012T ) and have been im- 
pressively used to mode l specific Be stars in c onsiderable 
detail (see, for exam ple, iCarciofi et al.l (|2009f ) for £ Tau). 
iCarciofi et al.l (|2012f ) have recently used the observed de- 
cline of the V-magnitude of 28 CMa during a disk dissi- 
pation phase to infer a value of a — 1.0 ± 0.2 for the disk 
viscosity parameter using detailed hydrodynamic and ra- 
diative transfer modelling. However this current paper is 
focused mainly on the geometrical origin of the correla- 
tion classes noted in the previous section and the influ- 
ence of the disk scale height and gravitational darkening 
of the central stars on these correlations. In addition, the 
models here are not meant to represent any specific Be 
star in detail. For these reasons, a very simple model was 
used to build the disk as a function of time (variations 
of this basic description are discussed at the end of this 
section) . The overall density of the disk was assumed to 
linearly increase over a time /3 to a maximum value of 
p max , namely 




The outer edge of the disk was assumed to expand lin- 
early with time with a fixed radial expansion speed v r 
as 

R d (t)=v r t. (5) 

For both spectral types, we took p max = 10~ 10 gcm~ 3 
and /3 = 365 days. For the B1V model, we took v r = 
5kms _1 and for the B5V model, we took v r = 3kms _1 . 
These choices ensure that the disk sizes expressed in 
terms of the stellar radius are the same for both mod- 
els. Hydrodynamical simulations sugge st outflow veloc- 
ities in visc ous disks of a few kms -1 (|Lee et al.l 119911 : 
lPorterl[l999l ). In addition, in the density law of Eq. ([1]), 
n was taken to be 3.5 and The was set to 0.6T c g. The 
growth of the disk with time for both models is given 
in Table [3] After one year, the disk around the B1V 
model reaches 1.05 AU in radius and has a mass of 
3.25 • 10- 9 M Q ; the B5V model reaches 0.63 AU with 
a mass of 6.50 ■ 10 -10 M Q . We note that both of these 
models result in quite substantial Be star disks, with 
low- inclination Ha equivalent widths approaching 40 A 
in emission. 



Table 2 

Disk density and radius as 
a function of time for the 
B1V and B5V models. 



Time 


PO 






(days) 


(gem' 


' 3 ) 


(R*) 


9.1 


2.510" 


12 


1.9 


18.2 


5.010" 


12 


2.8 


27.4 


7.510" 


12 


3.6 


36.5 


1.010" 


11 


4.5 


91.3 


2.510" 


11 


9.8 


127.7 


3.510" 


11 


13.4 


182.5 


5.010" 


11 


18.6 


237.3 


6.510" 


11 


23.9 


273.7 


7.510" 


11 


27.5 


310.3 


8.510" 


11 


31.0 


365.0 


1.010" 


10 


36.3 



Note. — The final disk 
masses after one year are 
3.25 • 10 -9 M (B1V) and 
6.50- 10" 10 M Q (B5V). 



To dissipate the disk as a function of time, a inner hole 
was evacuated in the final density model of the build 
phase at t = 365 days starting at the inner edge of the 
disk at 1 i?* . The outer edge of the hole was assumed to 
propagate outward at speed v r so that the disk sizes of 
Table [5] also represent the outer radii of the evacuated 
hole at time t + 365 days. Inside this hole, the den- 
sity was reduced by a factor of 10 3 . Thus the mass of 
the disk is reduced by a factor of 10 3 over the duration 
of the dissipation phase (one year) and the final density 
model for both disks was n = 3.5 and po = 10~ 13 g cm" 3 . 
This density is so low that the Ha line of the model 
is indistinguishable from the photospheric (absorption) 
line and the IR excess is < 1%. Thus the assumption 
is that the disk is lost "inside- out," a direction con- 
sistent with many observations (iRivinius "eTaLl l2l)01at 
IClark et al.ll2lMlMeTlland et al.ll2006lh The disk mass 
as a function of time for the B1V model is shown in 
Figure Q] Note that in the dissipation stage, the sharp 
density contrast at the outer edge of the evaluated hole 
was smoothed somewhat to improve the radiative trans- 
fer solution there. 

The calculatio n of the Ha line profiles was done with 
the BERAY code (|SigutJl2011f ). Here the radiative trans- 
fer equation is solved along a series of rays through the 
star-plus-disk system. For rays terminating on the stel- 
lar surface, phot ospheric LTE Ha profi l es wer e computed 
with the code of iBarklem fc Piskunovl (|2003l ). To obtain 
the V-magnitude of the system, a visual SED was com- 
puted fo r each model and the V-transmission curve of 
Johnson (196^) was used. 

Only a single bedisk model is required for the disk- 
building sequence as the disk can truncated at the ap- 
propriate outer radius in the BERAY code. However, an 
individual bedisk model is required for each time in the 
disk loss sequence. As the disk is evacuated inside-out, 
the growing size of the inner gap changes the photoion- 
izing radiation field reaching the remaining parts of the 
disk. Hence a consistent thermal structure for the disk is 
required in each case, and the temperature structure of 
the disk changes along the loss sequence. This is further 
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x 10"° 

3.5 1 




100 200 300 400 500 600 700 800 

Time (Days) 

Figure 1. Disk mass as a function of time for the B1V model. 
The maximum disk mass is 3.25 ■ 10 — 9 Mg at one year. The cir- 
cles represent the basic model described by Eqs. (J4j and J5j- The 
squares represent the constant po model and the diamonds, the 
flat n = 1 model, as described in the text. The disk mass as a 
function of time for the B5V model is very similar except that the 
maximum disk mass is 6.50 ■ 10 -10 Mg after one year. 

discussed in Section l4~2l 

In addition to the basic set of disk models described 
above by Eqs. Q and ([5]), two alternate assumptions 
about the growth of the disk mass during the build 
phase were considered. In the first variation, Eq. Q 
was replaced with the assumption of a constant base 
density, p (t) = pa, for all times, where pa was taken 
to be 10 -10 gcm -3 . In the second variation, the con- 
stant density p (t) — pa assumption was used along 
with the requirement that in the inner region of the 
disk where R < 5 i?* , the power-law index n in Eq. (fT|) 
was n = 1.0. An index of n = 3.5 was retained for 
R > 5 i?*. In this case, to give a model with a final disk 
mass consistent with the other two models, po was set to 
2.7- 10~ 12 gcm~ 3 . Figure [I] compares the disk mass as 
a function of time in the build phase for these two vari- 
ations with the basic model previously described. Both 
of these variations result in much more massive disks at 
early times but result in the same total disk mass at the 
end of the disk build phase at 365 days. 

4. RESULTS 

4.1. Basic Model 

To establish a reference set of models, the disks sur- 
rounding the stars were assumed to be in strict vertical 
hydrostatic equilibrium with no gravitational darkening 
of the central star. For typical stellar parameters, Eq. © 
predicts very thin disks with (H/R) ~ 0.05 in the inner 
regions. In addition, it is well known that Be stars are 
rapid rotators and hence gravitational darkening is po- 
tentially an important effect. This will be discussed in 
Section EU 

As discussed in the Section [3J stars showing a posi- 
tive correlation show a decrease in visual magnitude with 
an increase in the equivalent width (EW) of Ha. Stars 
showing an inverse correlation show an increase in visual 
magnitude with an increase in the EW of Ha. For the 
calculations in this report, instead of the visual magni- 
tude itself, the change in visual magnitude (AV) is used 
which is defined as the difference between the visual mag- 
nitude of a given model and the visual magnitude value 



for the last model of the disk loss phase at 730 days. As 
noted in the previous section, this last model shows no 
evidence of disk emission and can be taken to represent 
the isolated B star. 

Consider first the B1V and B5V models seen nearly 
face-on at i — 15°. Figure [2] shows the disk growth 
(1 — 365 days) and dissipation (366 — 730 days) phases 
of both models. Both spectral types behave in a very 
similar manner. The top panel of the figure shows the 
change in EW of Ha emission line as a function of timeJ3 
while the bottom panel shows the change in visual mag- 
nitude (AV). In the top panel, the Ha equivalent width 
increases with time as the density and the outer radius 
of the disk increase. It can be seen that the strength of 
Ha reaches its maximum, rs 30 A, by essentially ~ 100 
days, when the outer radius of disk has reached ~ 10 i?„ 
(Table [5]) . This is consistent with the observed sizes of 
Be star disks resolved with optical interferometry us - 
ing Ha (jTvcner et all 120051 iGrundstrom fc Giesl l200ll . 
When the disk dissipation stage starts at 365 days, there 
is a temporary increase in the Ha equivalent width for 
another ~ 15 days before entering a steep decline. This 
rather counter-intuitive result occurs because the equiv- 
alent width is a measure of the strength of Ha relative 
to the local continuum, i.e. 

EWS / (^) i ' 1 (6) 

where F v is the flux at frequency v and F c is the flux 
of the reference continuum. Because of the very large 
optical depths in Ha, it is formed over a wide region of 
the disk extending out to ~ 10 i?*. The optical con- 
tinuum, on the other hand, forms closer to the star. So 
in an inside-out dissipation scenario, the strength of the 
reference continuum is reduced more quickly than that 
in the line, and this leads to an increase in the equivalent 
width, even though F v declines during this phase. This 
'bump' is over by ~ 400 days, by which time the inner 
evacuated zone has expanded to ~ 3.5 i?*. 

The growth and dissipation phases do not mirror one 
another because of the inside-out nature of the dissipa- 
tion. In the growth phase, the base density (po) and 
outer radius (Rd) of the disk both grow linearly in time. 
However, during the dissipation phase, only the extent 
of the inner evacuated region grows linearly with time, 
leading to an inner cavity in the disk. By 730 days, the 
disk has p — 10 _13 gcm -3 and extends to 36 i?*. This 
disk is so rarefied that there is essentially no emission due 
to the disk in Ha, and its EW represents the (negative) 
absorption width of the photospheric line. Note that the 
first calculated model in the growth phase is after 9 days 
at which point the disk has grown to 1.9 i?* . This is 
the reason the initial equivalent width in the Figure [5] 
is slightly larger (i.e. affected by emission) than the last 
model in the dissipation sequence. 

The bottom panel of Figure [2] shows the change in 
the V-magnitude over the growth and dissipation phases. 
The reference model for the colour differences is chosen 
to be the last model in the dissipation sequence (730 
days) because, as noted above, it exhibits a pure photo- 

2 We adopt the non-standard convention that a positive EW 
denotes emission and a negative EW, absorption. 
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spheric spectrum which is uncontaminated by the disk. 
The sense of the visual magnitude changes are always 
brightening in the growth phase and dimming in the dis- 
sipation phase (with one brief exception noted below). 
As the system is viewed nearly face-on, there is little 
obscuration of the star by the disk, and hence the disk 
always adds to the brightness of the system. The very 
rapid decline in the V magnitude of the system at the 
onset of the dissipation phase reflects the above discus- 
sion about the formation of the optical continuum close 
to the star, within a few stellar radii. 

In the correlation classification scheme of Section [21 
both models (B1V and B5V) would be an instances of a 
positive correlation as the system brightens as the disk 
(and the Ha emission strength) grows. Note, however, 
there would be a brief window of an inverse correlation 
during the initial phase of the inside-out dissipation, as 
the system would grow fainter while Ha continued to 
increase during the 'bump' between 365 and 380 days. 

Another interesting way to look at this sequence of 
models is in a colour-colour diagram of (U — B) versus 
(B — V) where UBV are the familiar Johnson photomet- 
ric colours. The evolution of the entire sequence for the 
B1V and B5V models is shown in Figure [3] To calibrate 
the colours, the colours of the last model in the dissipa- 
tion phase (730 days) were forced to be the colours of 
a norm al B1V or B5V star as tabulated by iFitzgeraldl 
(1970). The additive constants were then applied to ev- 
ery model in the sequence. Note that as the first model 
is after 9 days of disk growth, it will not lie exactly on 
the main sequence reference point. For both modes, the 
growth track is from the main sequence, nearly horizontal 
in the figure, towards the supergiant sequence, resulting 
in models that have nearly constant (U — B) colours but 
progressively redder (B — V) colours. The tracks are such 
that the systems move from the V region, horizontally, 
towards the la region of the colour-colour diagram; this 
is the observed behaviour for positively correlated sys- 
tems in which movement in the colour-colour diagram is 
often categorized as the system "changing its luminosity 
class between V and la without any s ubstantial change 
of its spectral type" (jHarmancc 119831 ). This behaviour 
is a reflection of the nature of the additional continuum 
emission from the disk. Typically, the continuum excess 
for Be stars grows with longer wavelengths (because it is 
principally due to free- free emission). Hence the U colour 
is hardly affected, while the V colour is most affected (as 
noted above). As the star is essentially unobscured for 
i = 15°, this simple fact gives the horizontal motion in 
the colour-colour diagram. 

We now turn to the same set of models, B1V and B5V, 
but for a nearly edge-on disk seen at i — 89°. In this 
case, the disk can obscure the star's photospheric light, 
and this leads to new effects shown in Figure |U The 
top panel shows the change in the Ha EW, and it again 
grows with time as the disk is built. It starts at negative 
values this time (i.e. net absorption) due to the much 
weaker disk emission at i = 89 degrees and the effect 
of shell absorption. The photospheric equivalent widths 
are represented by the equivalent widths at 730 days. In 
this series, the Ha equivalent width grows steadily with 
time and reaches a maximum in emission of ks 6 A at 
365 days. In the edge-on case, the maximum Ha EW 
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Figure 2. Ha equivalent width (A, top panel) and change in 
V-magnitude (bottom panel) as a function of time for the B1V 
(solid lines) and B5V (dashed lines) models seen nearly pole-on at 
i = 15°. All models assumed a thin, hydrostatic disk with H set 
by Eq. and gravitational darkening of the central B star was 
not included. The disk grows for one year and then dissipates over 
one year as described in the text. 

is much smaller than in the face-on case, as expected 
from the much smaller projected area of the optically 
thick disk on the sky. The bottom panel shows the ef- 
fect on the V magnitude. For times less than about 100 
days, the system is fainter than the star itself, reflecting 
the obscuration by the disk. However, the effect is very 
small (a few hundredths of a magnitude), and this is a 
direct consequence of the very thin nature of a purely 
hydrostatic disk. 

In the colour-colour diagram corresponding to the i = 
89° models (Figure |3|), more inclined tracks are seen. In 
this case, the U colour is affected by absorption by the 
disk, making it larger and hence (U — B) smaller. Obser- 
vationally, inversely correlated systems "move along the 
main s equence changin g its spectral but not luminosity 
class" (|Harmaned [T983) . This behaviour is clearly repro- 
duced for the B5V model at i = 89° during its disk build 
phase. However, the B1V model does not move along 
the main sequence, and this may be related to the very 
weak inverse correlations predicted by this set of basic 
models with hydrostatic disks. 

4.2. Enhanced Disk Scale Heights 

In the previous section, we presented the optical signa- 
ture of a disk build and dissipation sequence. The case 
of positive correlation is well represented by the i = 15° 
models of both spectral types. The overall V magnitude 
of the system increases by approximately half a magni- 
tude over the build phase of the disk. In addition, the 
movement in the colour-colour diagram is as observed: 
towards the giant sequence at nearly constant spectral 
type. 

However, the case of an inverse correlation is less well 
represented. In the edge-on model, for which an inverse 
correlation is expected, only a very small magnitude in- 
crease is seen during the initial build phase of the disk, 
amounting to a few hundredths of a magnitude. Ob- 
served inversely correlated systems typically show mag- 
nitude changes of a tenth of a magnitude or more. Be- 
haviour in the colour-colour diagram is more as observed; 
there is a tendency to evolve along the main sequence, 
changing the system's spectral type. 
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Figure 3. The movement of the B1V (solid lines) and B5V 
(dashed lines) models in an optical colour-colour diagram as the 
disk grows (filled blue circles) and dissipates (open red circles). The 
models assumed a thin, hydrostatic disk and gravitational darken- 
ing was not included. The models are viewed nearly pole-on at 
i = 15°. The circles grow and diminish in size following p a (t)- 
The squares are the observed normal (star only) colours for m ain 
sequence (V) and supcrgiant (la) stars from Fitzgerald (1970). 




Figure 4. Same as Figure[2]but for edge-on disks seen at 



solid: B1V i = 89° 
dashed: B5V i = 89° 
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Figure 5. Same as Figure[2]but for edge-on disks seen at ■ 



One feature of the reference models leading to the very 
small magnitude increase during edge-on observation is 
the very thin nature of the disks. Almost all calcula- 
tions of Be star disks, either hydrodynamical or ther- 
mal/radiative, assume a geometrically thin disk with the 
scale height set by vertical hydrostatic equilibrium. As 
previously noted, this scale height is the ratio of the lo- 
cal sound speed to the local orbital speed. As Be star 
disks are relatively cool and the orbital speeds are Keple- 
rian and large, the disk scale heights are predicted to the 
very small; near the star, the scale height is often only a 
few percent of the stellar radius. Some have questioned 
whether this assumption is really correct, and there is 
some circumstanti al evidence to support enhanced inne r 
disk scale heights (jArias et al.ll2006l ; IZorec et al.ll2007bD . 

To examine the sensitivity of the predicted magnitude 
changes to the inner disk scale height, we have adopted 
a simple model in which the disk scale height is not al- 
lowed to fall below a minimum value, (H / R) m i n , which 
we take to be 0.2. The effect on the disk scale height of 
the B1V model is shown in Figure [6] As can be seen, 
this assumption results in a disk of nearly constant scale 
height (H/R) = 0.2 until « 30 stellar radii. Near the 
star, the implied scale height enhancement is a factor of 
ps 5. Given this modification, the temperature structure 
of the disk, the Ha EWs, and the optical UBV magni- 
tudes and colours were recomputed. 

One subtlety of the enhanced scale height models is 
that simply increasing the scale height over the refer- 
ence model will result in disks that are more massive. 
To avoid this effect, we have reduced the po of the en- 
hanced models at each time in Table © by the amount 
required to match the disk mass of the reference model. 
Hence the disk mass of the enhanced scale height models 
as a function of time match those given in Figure [TJ We 
note that the total emission measures of the disks now 
agree much better as well (within 30%). Another impor- 
tant effect of this reduction in p is that the temperature 
structure of the enhanced scale height models are now 
much closer to those of the basic reference models. This 
is shown in Figure [7] Without the reduction in p , the 
enhanced scale height model for B1V has a much more 
extensive inner cool zone in the equatorial plane as com- 
pared to the reference hydrostatic model. The enhanced 
scale height model with reduced po shows much better 
agreement. 

In this initial set of enhanced scale height calculations, 
we ignored gravitational darkening; this is potentially an 
important effect and it will be modelled in Section I4.4I 
Hence the results here can be used for Be stars rotating 
at less than » 70% of their critical velocity. Instead of 
presenting the entire sequence of models, as before for the 
reference models, we will focus solely on the correlation 
between the optical magnitude and Ha equivalent width 
during the disk build sequence from I to 365 days. 

The results are shown in Figure |S] for the B1V mod- 
els and Figure |9] for the B5V models. The contrast with 
the reference hydrostatic models is striking. The thin, 
hydrostatic models show only a very small decrease in 
the brightness of the system while the enhanced scale 
height models clearly show strong inverse correlations 
for inclinations above about « 75°. Perhaps the most 
significant difference is that increases in the visual mag- 
nitude of the system occur even for the largest values of 
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Figure 6. Top panel: hydrostatic scale height (Eq. [3J for the 
B1V model assuming T HE = 15 000 K. The scale height (H/R) at 
1 H* is predicted to be ~ 0.04. (H / R) m \ n = 0.2. adopted for some 
models in this work, is shown as the dashed line. Bottom panel: the 
disk scale height enhancement over the hydrostatic value obtained 
adopting a minimum value for (H/R) of 0.2. The dotted line is 
only a guide for the eye. 




Figure 7. Comparison of the equatorial disk temperature between 
the reference B1V model (solid line) and a model with (H j R) m i n 
set to 0.2 (dashed line). Also shown is the temperature structure 
of a (H I -R) m i n = 0.2 model with po reduced so that the total disk 
mass equals that of the reference B1V model (dash-dot line). 

Ha EW. In the pure hydrostatic models, all of the curves 
turn over at intermediate EW, and the system becomes 
brighter. From the enhanced scale height models, the 
dividing viewing inclination between positive and nega- 
tive or inverse correlation is seen to be k 75°. Hence 
the choice of (H/R) min = 0.2 seems to give an i crit in 
agreement with the rough estimate of Section [2] based on 
the observed numbers of positive and negative correla- 
tions. Nevertheless, this should not be over interpreted; 
a more homogeneous sample of correlations coupled with 
detailed modelling of individual stars is required for a 
firm conclusion on enhanced disk scale heights for Be 
stars. 

4.3. Alternate Disk Growth Models 

One feature of the disk growth model implied by 
Eqs. (01 and §5§ is that the disk density during the early 
growth phases, t < 100 days, is very low, and hence there 
is not a large impact on the system's spectral energy dis- 
tribution. To see the effect of a much denser disk during 
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Figure 8. Relation between the change in V-magnitude and the 
Ha equivalent during the disk building phase (1-365 days) for the 
B1V model. Each curve is labelled by the inclination of the system 
to the line of sight. The left panel shows the case of pure hy- 
drostatic disks and the right panel shows the case of models with 
enhanced vertical scale heights. 
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Figure 9. Same as Figure l8l except for the B5V model. 

the early growth stages, we return to the two alternate 
disk growth assumptions outlined at the end of Section|3l 
Note that both of these alternate models assume a thin 
disk with the vertical scale height set by hydrostatic equi- 
librium via Eq. ([2]); these are not enhanced scale height 
models. 

In the first alternate model with the base disk density 
set to a constant p Q (t) = 10~ 10 gcm~ 3 , the outer ra- 
dius of the disk simply grows with time. For this model, 
we have recalculated the change in V-magnitude and 
Ha equivalent width during the growth phase (t < 365 
days), and the results are shown in the left-hand panel 
of Figure [10] for the same viewing angles as the previ- 
ous figures. As can be seen from the figure, the A V at 
early times is much larger because of the denser disk. 
The period of inverse correlation, however, is confined 
to the very initial times when the Ha EW is still small. 
For later times, and larger values of the Ha EW, the 
correlation flattens out and then turns over for all view- 
ing angles. The case of the B5V model is shown in the 
left-hand panel of Figure [TT] and is very similar to the 
B1V model, except that changes in the V magnitude are 
somewhat smaller. 

In the second alternate disk growth model, a constant 
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Figure 10. Same as Figurel8lfor the B1V mod el bu t implementing 
the alternate disk growth scenarios of Section 14.31 The left panel 
is the constant po disk, while the right panel has constant po and 
n = 1 in the inner region. 
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Figure 11. Same as Figure ITOl but for the B5V model. 

base disk density was again used, and in addition, a 
much shallower density drop-off in the equatorial plane 
for Rd < 5 i?* was assumed by setting n = 1 in Eq.[T]over 
this range. To ensure that the fully built disk had the 
same mass as the previous models, the base disk density 
was taken to be a constant p (t) = 2.7 • 10~ 12 gcm~ 3 ; 
the disk mass as a function of time is also shown in Fig- 
ure [TJ The results in the A V-Ha diagram are shown 
in the right-hand panel of Figure [TU] for the B1V model. 
In this case, weak inverse correlations are predicted for 
inclinations i > 80°, but the dimming of the system is 
small, w 0.05 mag. In addition, like the previous 
case of p (t) — 10~ 10 gcm -3 , the inverse correlation is 
confined to small Ha EW; the relation flattens out and 
turns over into a positive correlation for larger EWs. A 
very similar result is obtained for the B5V model, as 
shown in the right-hand panel of Figure 1111 

4.4. Gravitational Darkening 
Be stars a re known to be rapid rot ators (IPorterl 1996; 



lYudinl 120011 ). as first suggested by iStruvd ([193 If ), al 
though how close they are as a population to crit- 
ical rota t ion i s still uncertain dTownsend et al.l l2004t 
iCranmerl l2005t iMeilland etldl |2012(). Rapid rotation 
leads to gravitational darkening of the stellar surface 
in which the temperature varies with latitude, and the 



shape of the star is n o longer spherica l. According to 
von Zepiel's theorem ()von Zeipell 119241 ) . the local effec- 
tive temperature should be proportional to the effec- 
tive gravitational acceleration at each latitude, and this 
effect causes the equatorial regions of a rapidly rotat- 
ing star to be much cooler than the poles. Gravita- 
tional darkening, both as the distortion of the stellar 
surface and as a variation of temperature over the stel- 
lar surfac e, now has amp le, direct interferometric ver- 
ification (jvan Belle! I2012D . These studies have lead to 
the conclusion that while gravity darkening is present, 
von Zepiel's th eorem as applied to stellar atmospheres by 
iCollinsI (119661) overestimates the temperat ure variation. 
Recently. lEspinosa Lara fc Rieutordl (|2011l ) have derived 
a new version of gravitational darkening that seems in 
better agreement with observations; we adopt their for- 
malism. 

Gravitational darkening is potentially an important ef- 
fect in the current work, particularly for the case of an 
inverse correlation because the blocked equatorial region 
is now much cooler and hence fainter. One might expect 
that the size of the inverse correlation will depend on the 
extent of gravitational darkening. 

We have used the formalism of 

lEspinosa Lara fc Rieutordl (|2011[ ) to compute the 
variation of temperature with latitude over the stellar 
surface. Given the exploratory nature of the present cal- 
culations, however, we have neglected the distortion of 
the stellar surface. A useful term for this approximation 
is spherical gravitational darkening (SGD). We have also 
neglected the effect of gravitational darkening on the 
thermal structure of the circumste llar disk. Thi s latter 
effect is extensively discussed by iMcGill et al.l (|201 ll ) 
and even at critical rotation, the effect would be at most 
a 15% reduction in the average disk temperature. 

Figure [12] shows the effect of SGD, with the central 
star rotating at 95% of its critical velocity, on the B1V 
model, both with a thin, hydrostatic disk (left panel) 
and the thicker (H/R) m \ n = 0.2 disk (right panel). In 
the case of the thin, hydrostatic disk, almost all trace of 
the inverse correlation is gone. In the (H/R) m - m = 0.2 
case, the inverse correlation is now seen only for higher 
inclinations (i > 80°) and only for latter times in the 
disk-building process. A very similar result is seen in 
Figure [H for the B5V model.' 

As a result, the existence of an inverse correlation de- 
pends strongly on the extent of gravitational darkening 
and hence, on the rotation rate of the central B star of 
the system. Note that these results would only contra- 
dict the existence of inverse correlations if all Be stars 
were rapid rotators. While the exact rotational distribu- 
tion of Be stars remains unclear, it seems unlikely that 
all of them are critical rotato rs, at least for the earlier 
spectral types (|Cranmeril2005D . 

5. DISCUSSION 

In this work, simple parametrized models for Be star 
disk growth and dissipation are investigated to see if they 
can reproduce the known classe s of positive and inverse 
correlations discovered by Harmancc (1983). We find 
that these simple models can demonstrate both positive 
and inverse correlations between the system's change in 
visual magnitude and the magnitude of the Ha emission. 
We also demonstrate the scale height of inner circum- 
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Figure 12. Same as FigurelBlfor the B1V model but implementing 
spherical gravitational darkening assuming the star is rotating at 
95% of its critical velocity. 
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Figure 13. Same as Figure H"2l but for the B5V model. 

stellar disk and the extent of gravitational darkening of 
the central star's surface play key roles in controlling the 
magnitudes of the predicted correlations, particularly the 
inverse ones. 

What do these calculations demonstrate? It is impor- 
tant to keep in mind that only parametrized, "toy" mod- 
els of disk growth and dissipation were used, and hence it 
is difficult to draw general conclusions. However the cal- 
culations performed suggest that inverse correlations are 
more easily predicted for Be star disks with enhanced 
scale heights around B stars that are not significantly 
gravitationally darkened. Because the predicted correla- 
tions depend so strongly on these two aspects, gravita- 
tional darkening and disk scale height, it would be highly 
desirable to analyze a well determine set of contempora- 
neous magnitude and Ha measurements through either 
a Be star growth or dissipation phase to determine to 
what extend the disk scale height and/or extent of grav- 
itational darkening can be constrained. In the analysis 
of real observations, a much more realistic hydrodynam- 
ical description of disk growth wou ld be preferable and 
such models are now available (e.g. ICarciofi et al.ll2012l : 
Haubois et alllMl) . 
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